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Abstract 

The random process of the tree evolution with continuous time 
parameter has been investigated. By introducing the notions of liv- 
ing and dead nodes a new model of random tree evolution has been 
constructed. It is assumed that at t = the tree is consisting of a 
single living node (root) from which the evolution may begin. This 
initial state is denoted by Sq. At a certain time instant r > the 
root produces v > living nodes connected by lines to the root which 
becomes dead at the moment of the offspring production, r and v 
are random variables with known distribution functions. In the evolu- 
tion process each of the new living nodes independently of the others 
evolves further like a root. The state St of the tree at time instant 
t > is defined by the actual numbers of living and dead nodes. It 
has been shown that the generating function of the probability to find 
a tree in a given state St at t > 0, assuming that at t = it was in the 
state So, satisfies a non-linear integral equation. Analyzing the time 
dependence of the expectation value of the number of living nodes it 
has been recognized that three completely different types of evolution 
exist depending on the average number q\ of nodes produced by one 
living node. If q\ < 1, then the evolution is subcritical, if q\ = 1, then 
it is critical, while in the case of q\ > 1 is supercritical. It has been 
proved that in the case of subcritical evolution there is a time point 
at which the variance of the number of living nodes is maximal. By 
using the the generating function of the probability to find at time 
moment t > a given number of nodes in the tree, the time depen- 
dence of the expectation value and the variance of tree size has been 
calculated. The joint distribution function of the numbers of living 
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and dead nodes has been also derived, and the correlation between 
these node numbers has been determined as a function of time pa- 
rameter. It has been proved that the correlation function converges 
to V3/2 independently of the distributions of v and r when qx — ► 1 
and t — > 00. 

PACS: 02.50.-r, 02.50.Ey, 05.40.-a 

1 Introduction 

There is a growing interest in randomly evolving graphs since a large number 
of practically important problems can be analyzed in this way. In the last 
few years so many papers have been published in this field that any list of 
references would be far from complete by any standards. It is very fortunate 
that recently two outstanding review papers 0, [@] have been published and 
thus the author of the present paper does not feel obliged to cite the large 
amount of important references. 

There are many different models for the random graphs and it would be 
difficult to give a precise classification of them.Q Probably the oldest model 
has been proposed by Erdos and Renyi ||, [Q, ||. In this model the initial 
state of the graph consists of fixed number n of nodes that will be randomly 
connected during the evolution. The time is discrete and counted by the 
number of successively appearing lines. At time t — 1 one of the f ™) possible 
lines defined by the nodes is chosen out assuming that each of these lines has 
the same probability to be chosen. At time t = 2 the second line is chosen 
out from the (™) — 1 possible lines all these being equiprobable. Continuing 
this process at time t = N a random graph is formed which is consisting of 
n fixed nodes and N lines selected randomly. If n is a large positive integer 
converging to infinity and N is increasing from 1 to \ Z\ then the evolution 
of the graph passes through several (say five) clearly distinguishable phases. 

One has to mention that random networks with fixed number of nodes 
usually show the so-called small-world effect, i.e. the average shortest path 
is small. This type of graphs can be created by randomly rewriting some 
of the lines of a regular ring graph. Watts and Strogatz || noticed many 

1 Throughout this paper we will use the terms "node" and "line" (sometimes "link") 
for the terms "vertex" and "edge", respectively. The degree of a node is defined by the 
total number of its links, and the out-degree of a node is nothing but the number of its 
outgoing lines. 
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interesting properties of these graphs called small- world networks and a large 
number of papers QJ, [Q was published recently on this subject. 

There is another - more realistic - model |J for the random evolution of 
graphs. In this model the number of nodes is not fixed, but to the contrary 
of the Erdos and Renyi model it is increasing in each discrete, equidistant 
time moment by addition of a new node connected randomly to one or sev- 
eral of the nodes already existing in the graph. The probability of a given 
connection between the new and one of the old nodes may depend on the 
number of links due to the node existing already. Many variants ||10H , [11 



of this scale-free model have been recently investigated by using different in- 
ventive approximations. The one of the most intensively studied problem is 
the chracter of the degree distribution of nodes in graphs evolving randomly 
with various preferential linking. The number of exactly solvable (or solved) 
problems, however, is not very much, therefore, it seems to be worthwhile to 
study models which can be treated exactly. 

It is known for long time that the random evolution of a tree | corresponds 
to a Galton- Watson process |12| which describes the formation of a popula- 
tion which starts at t — with one individual and in which at discrete time 
moments t = 1,2,... each individual independently of the others produces 
a random number of offspring with the same distribution. There are many 
interesting papers on this subject but it is hardly to find articles about 
the random processes of tree evolution with continuous time parameter. The 
purpose of is this paper is the study of these special random trees with the 
hope that it brings about new insight into the nature of branching processes. 

This paper is organized as follows. In Section 2 the exact description of 
the random process is given, while in Section 3 the derivations of the basic 
equations for the generating functions can be found. The properties of aver- 
age characteristics (expectation values, variances and correlation functions) 
are discussed in Section 4, and, finally, the conclusions are summarized in 
Section 5. 



2 Description of the random process 

We will study the random evolution of such trees which are consisting of 
living and dead nodes as well as lines connecting the nodes. The living node 

2 A randomly evolving tree is a connected graph containing no cycles, and growing from 
a single node (root) according to well-defined rules. 



3 



is unstable, and after a certain time r may become dead. The lifetime r is 
random variable, and denote by V{r £ (t,t + dt)} = dT(t) the probability 
of finding r in the time interval (t, t + dt). The dead node is inactive, cannot 
change its state, and is unable to interact with the other nodes. 

The evolution process can be described as follows. Let us suppose that 
at time instance t — the tree consists of a single living node called root 
and denote by <S this state of the tree. At the end of its life the root 
creates v = 0, 1, 2, . . . new living nodes and after that it becomes immediately 
dead. The number of the offspring v is random variable, and the probabilities 
V{v — k} — fk, k = 0, 1, 2, . . . are supposed to be known. The new living 
nodes are promptly connected to the dead one and independently of the 
others each of them can evolve further like a root. It means that the tree 
from its initial state So may evolve further by two mutually exclusive ways: 

• the root does not die, and thus the tree remains in its initial state 
during the whole time interval (0,t), or 

• the root does die in an infinitesimally small subinterval (t', t' + dt') e 
(0, t) and produces j > new living nodes. Each of these nodes will 
evolve further in the remained time interval (t',t) independently of the 
others like the root. 

Fig. [I] shows the initial state of a tree consisting of a single root and 
illustrates four of the possible states which may be produced by the dying 
root. The probabilities of producing 0, 1, 2, 3, . . . living nodes are denoted by 

fo, fl, f'2, f3, 

Fig. Q illustrates one of the realizations of a random tree at the moment 
t > 0). One can see in Fig. |2] three white circles denoting the living nodes 
which are capable for the further evolution and a large number of the black 
circles denoting the dead nodes which are unable to produce new living nodes. 
It is to mention that each of the offspring is connected by a line to the dead 
precursor. If a dead node has k outgoing lines then it will be called node of 
kth out- degree. It is obvious that a living tree must have living nodes, while 
a dead tree consists of dead nodes only. 

In order to simplify our consideration we will assume that the distribution 
function of the lifetime r of living nodes is exponential, i.e. 

dT(t) = e- Qt Q dt, (1) 
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Figure 1: First step of the tree evolution. The new state arises from the initial state 
which consists of a single root. Four of the possible new states are illustrated in Fig. The 
living nodes are denoted by white and the dead ones by black circles, fo, fi, fi, fs, ■ ■ ■ are 
the probabilities of producing 0, 1, 2, 3, . . . nodes by the dying root. 

where 1/Q = t is the expectation value of the lifetime. In this case the 
evolution becomes a Markov process. It is assumed that the probabilities 



are the same for all living nodes. By introducing the generating function 



V{v = k} = f k , fc = 0,l,2,... 



(2) 



oc 



?0) = E /* 



(3) 



fc=0 



we can define the factorial moments of v in the following way: 




J 2 = 1 




For the expectation value and the variance of v we have 

E{*4 = qi and D 2 {z/} = q 2 + qi - q\. 



(4) 



For the sake of later use we cite the relation 




(5) 
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Figure 2: One of the realizations of a random tree at the time moment t > 0. The nodes 
denoted by white circle are capable for the further evolution, while those denoted by black 
could not further evolve. The horizontal dashed lines indicate the time instances at which 
nodes were produced. 

In order to characterize the random process of tree evolution we introduce 
the random functions (J,e(t) and fi d (t) which are the actual numbers of the 
living and dead nodes, respectively, at time instant t > 0. It seems to be 
worthwhile to define the random function = ^(t) + fj, d (t) which gives the 
total number of nodes, i.e. the actual value of the tree size at t > 0. Finally, 
let us denote by fJ, e (t) the number of lines in a random tree at time instant 
t > 0. It is obvious that the equality (J> e (t) — /i(t) — 1 is true with probability 
1 for all t > 0. 

In the next Section we shall derive the basic equations for the gener- 
ating functions of probabilities of the following random events: {fJ>e(t) = 
n £ }, {/i d (t) = n d }, {ne(t) = n e , /i d (t) = n d } and {/i(t) = n}, provided that 
at t — the tree was in the state «S . 

3 Derivation of the basic equations 
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3.1 Generating function equations 

3.1.1 Probability distribution of the number of living nodes 

The living nodes capable to produce new nodes at the end of their life are 
responsible for the evolution of the tree. It seems to be interesting to have 
an insight into the random process of the formation of living nodes. Let us 
denote by 

V{fi e (t)=n\S }=p( e \t,n) (6) 

the probability that the number of living nodes fie(t) is equal to n at the 
time instant t > provided that at t — the tree was in the state Sq. By 
exploiting the rules of the random evolution described in Section 2 it can be 
easily shown that 



-Qt 



S n ,i+ 



Q 



-Q(t-t') 



,71 j 



k=l niH \-n k =n j=l 

and by using the probability generating function 



dt' 



g^(t,z) = Y,P ie) (t,n) 



n=0 



we can prove that g^ e '(t, z) satisfies the integral equation 



9 {i) (t,z) 



-Qt 



z + Q 



e -Q(t-f) q 



g {e) (t',z)] dt'. 



(7) 



It is an elementary task to prove that the integral equation ([?[) is equivalent 
to the differential equation 



dt 

with the initial condition 



Qg^(t,z)+Qq gW(t,z 



(8) 



limc/ w (t, z) 



z. 



From the equation (0) or (H) the time dependence of the expectation value 
and the variance of the random function /^(t) can be easily calculated. 
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3.1.2 Probability distribution of the number of dead nodes 

A living node producing offspring becomes dead, and we are interested in the 
probability that the number of dead nodes Hd{t) is equal to n at time instant 
t > provided that at t — the tree was in the state So- By introducing the 
notation 

V{fi d (t)=n\S }=p (d \t,n), (9) 
we can easily derive the following equation: 



P 



(t, n) = e~ Qt 5 nfi + 



Q I t 



-Q(t-t') 



/ <w + £/* E I[p (d) (*>;) 

k=l niH h«fc=n— 1 j=l 



and by using the probability generating function 



g^ d \t,z) = Y,P {d \t 



id), 



n) z 



n=0 



we can prove that g^ d \t, z) satisfies the integral equation 



dt', 



(10) 



g id \t,z) =e~ Qt + zQ 



Qt 



-Q(t-f) 



g {d \t',z)] dt'. 



It is obvious that the equation ( PH| ) is equivalent to the differential equation 

dgW(t,z) 



dt 

with the initial condition 



QgW(t,z)+zQq g^ d \t,z 



]imgW(t. z) = 1. 



(12) 



3.1.3 Probability distribution of the total number of nodes 

Denote by /i(t) = fJ,e(t) + ^dif) the sum of the numbers of living and dead 
nodes at t > 0. Let us determine the probability of the random event = 
n} provided that at t — the tree was in the state Sq. Introducing the 
notation 

V{fj.(t) =n\S } =p(t,n) (13) 
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and denning the generating function 

oo 

g(t,z) = y Ep(t,n)z n , (14) 

n=0 

it can be easily shown that g(t, z) satisfies the integral equation 

g(t, z) = e~ Qt z + zQ t e ~ Q ^ q \g(t', z)\ dt' . (15) 

Jo 

For the sake of completeness we are deriving from Eq. fll5|) the equivalent 
differential equation: 

^ = ~Q 9(t, z) + Qzq [git, z)] 
the initial condition of which is 

lim q(t, z) = z. 

no v ' 

3.1.4 Joint distribution of the numbers of living and dead nodes 

We shall define now the probability 

V{fi e (t) =ne,n d (t) = n d \S } = p {e > d) (t,n e ,n d ), (16) 

where ri£ and n d are a non-negative integers. It is clear that p^' d \t,n£,n d ) 
gives the probability that at the time instant t > the number of living 
and that of dead nodes in a randomly evolving tree are equal to ne and n d , 
respectively, provided that at t = the tree was in the state <So- By exploiting 
the rules of the random evolution described in Section 2 it can be shown that 

p^ d \t,n e ,n d ) = e~ Qt S neA S ndfi + 
Q J o e-^j/o^o^i + ^/fcy^^^n^J dt, 

where 

V^ d \t>, n t , n d ) = Yl Eft P M V, ntJ, 

n«,lH \-ne,k= n e ™d,H \-n d! k=n d —l j=l 
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Introducing the probability generating function 



oo oo 



g^ d \t,z e ,z d ) = J2 E^M.ni) 

n e =0 n d =0 



Zo Z. 



d > 



(17) 



and taking into account the generating function defined by Eq. (^) we obtain 
the fundamental equation in the form: 



g^ d \t } z e ,z d ) = e-^ z e + z d Q 



Qt 



e -Q(t-f) q 



g^(t',z £ ,z d ) dt', (18) 



which can be rewritten as a nonlinear differential equation 



d \t,z t ,z d ) 
dt 



-Q g^ d \t, z t , z d ) + z d Qq gM(t, z t , z d ) (19) 



with initial condition 



\img {l4 \t,zz,z d ) = z t . 



One has to note that the former equations (|7]) and (ITTp can be immediately 
derived from the integral equation (Pf). 



3.1.5 Joint distribution of the numbers of living nodes and lines 

As it has been already mentioned the living nodes are responsible for the 
evolution of the tree. It seems to be interesting to have a deeper insight into 
the interplay between the living nodes and lines during the tree evolution. For 
this purpose let us determine the probability that the number of living nodes 
Ht{t) born in the time interval (0, t) and the number of lines fi e (t) produced in 
the same time interval are equal to and n e , respectively, provided that at 
the time instant t — the tree was in the state Sq. Introducing the notation 



V{fj,i(t) =ne,n e (t) = n e \S } = p {e ' e \t, ne,n e ), 
we can easily derive the backward equation 

p {e ' e) (t,n e ,n e ) = e~ Qt 5 ni , x 5„ e , + 



(20) 



Q I t 

10 



-Q(t-t') 



fo <W,o K e ,G + ^2fk M k (t', n e , n e ) 
k=i 



dt' 
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where 

M k (t',ni,n e ) = 

The probability generating function 



oo oo 

g<M(t, zi, z e ) = £ E P ( *' e) (t, nt, n e ) zT z, 

ri£=0 n e =0 



E E U P^raj^j). 

iH Yn l>k =n t n e ,iH \-n e ^ k =n e -k j=l 



e 5 



(21) 



we obtain 



g^(t, z e , z e ) = e-^ z i + Q f e"«<«-*> q[z e g^\t', z e , z e )\ dt' . (22) 

Jo 

This is an important equation from which one can easily derive relations for 
the time dependence of factorial moments fie(t), and // e (t), and calculate the 
correlation function between them. The particular variant of the Eq. (p2"D, 
namely: 



g^\t,z e = z,z e = l)=g^(t,z) = e-^z + Q I e 



-Qt 



-Q(t-t') 



9 {i) (t',z) 



dt' 



is equivalent to Eq. ([7|), and by substitutions zi — 1 and z e = z one can 
obtain the integral equation 



g^\t,z e = l,z e = z)=g^(t,z)=e- Qt + Q I e 



(e)( 



-Qt 



,-Q(t-t) 



zg( e \t',z)] dt' 



determining the generating function g^ e \t,z) = E{z Me }. 

The differential equation corresponding to (E21) has the form: 



dg^(t,z t ,z e ) 
dt 

with the initial condition 



Qg^ e \t } z e ,z e ) + Qq z e gW(t,z t ,z e ) 



limg^ e) (t, 



Z e ) = Z £ . 



11 



3.1.6 Remarks on generating function equations 

In this Section we derived integral equations for the following generating 
functions: 

gW(t,z) = E{z' H ®}, g (d \t,z) =E g(t,z) =E{z"®} 



and 



g^(t, z e , z d ) = E {tf® #®} , g^(t, z e , z e ) = E [z^ #<«>} . 

For the sake of simplicity let us denote by j(t, z) f\ any of these generating 
functions. It can be shown |L4] that the integral equation for j(t, z) has just 



one solution satisfying the conditions 

|7(t,^)|<l, V \z\<l and lim ^(t,z) = l. 

zfl 

It is important to mention that if lim 2 -f 1 7(£, z) < 1, then the number of 
the tree elements (nodes, lines) increases so rapidly that it reach 00 in a 
finite time interval. In our further considerations, however, this explosion- 
like evolution will be excluded. 



4 Average characteristics 

4.1 Expectation values and variances 

One of the simplest way to study the characteristic properties of tree evolu- 
tion is to analyze the time dependence of factorial moments of the following 
random functions: 

He{t), fJ>d(t), and n(t). 

For the sake of simplicity let us denote by x(t) any of these random functions 
and define the generating function 

7 (M)=E{3*W}. (23) 
According to the Abel's theorem the jth factorial moment 

>7(t, z 



E{x(*) \x(t) -!]••• [X(t)-J + 1]} = mjit) 



dzi 



, (24) 

2=1 



3 The variable z may have two components. 
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exists if and only if the series 

oo 

J2V{ X {t) = n\S }n j , j = 1,2,..., 

n=0 

converges uniformly at any finite time instant t > 0. It is well-known that 
the simple moments 

E{[ X (tW} = m^(t) 
can be expressed by the factorial moments ]15] by using the relation 

m { *\t) = Y,S{j,i)m i {t), j = l,2,... , (25) 

i=0 

where S(j,i), i = 0, 1, . . . , j are the Stirling numbers of the second kind. 

In the further considerations we are going to deal with two important aver- 
age characteristics of the random function x(t), namely with the expectation 
value E{x(t)} = mi{t) and the variance D 2 {x(t)} = m2(t)+mi(t) [1— mi(t)). 

In order to calculate these average characteristics we do not have to know 
the distribution function of z/, it is enough to know the values of the first and 
second factorial moments of v only. In this case the distribution of v we call 
arbitrary (a). However, if we are interested in solving the generating function 
equations, then we should have the exact form of the distribution of v. In 
the sequel, we will use two simple distributions, namely the Poisson (p) 

V{v = k} = |, 

and the geometric (g) 

distributions. In the case of the Poisson distribution we have 

E{*4 = qi and D 2 {z/} = q h 
while in that of the geometric one 

E{u} = qi and B 2 {u} = qi {l + qi ). 
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4.1.1 Living nodes 

Let as investigate the expectation value and the variance of the number of 
living nodes at the time instant t > 0. The equation for the expectation 
value E{/^(t)} = m!p(t) can be immediately derived from the Eq. We 
obtain 

mf{t) = e~ Qt + Qq x f e- Q ^ mf\t') alt', 

Jo 

which is a typical renewal equation. After elementary calculations the solu- 
tion is obtained in the form: 

mf\t) = e~^- qi)Qt , (26) 

from which we see that the average number of the living nodes has the fol- 
lowing asymptotic properties: 

0, if qi < 1, 
}im = { 1, if gi = 1, (27) 
oo, if qx > 1. 



t^oo 



The evolution of the tree will be called subcritical, if gi < 1, critical, if q\ = 1, 
and supercritical, if q\ > 1. 

In order to calculate the variance D 2 {/i^(t)} we need the second factorial 
moment nvp (t) which can obtained by solving the integral equation 



o 



mf{t) = Q I e^C*-*') \ q t mf{t') + q 2 \mf\t') 
It is easy to show that 



alt'. 



( [1 - e-^-i^] e -^)\ if qi ^ l, 
mf{t) = \ ^ (28) 
{ Q2 Qt, if qi = 1, 

and by using the relation 

D 2 {M*)} = mf{t) + mf\t) [1 - mf\t)} 

we have 

Hkzll!} [i _ e -d-9i)Qt] e -(i- ffl )Ot j if gi ^ i, 
D 2 {/Mf)}-<1 91 (29) 
92 Q*, if <?i = 1, 
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where E{(z/ — l) 2 } = (1 — g x ) 2 + ~D 2 {v}. In the case of subcritical (qi < 1) 
evolution the variance has a maximum at 




20 40 60 80 100 

time in Q" 1 units 



Figure 3: The variance of the number of living nodes vs. time parameter Qt in subcritical 
(qi = 0.95) evolution. The curves g, p, and a correspond to the geometric, Poisson, and 
the arbitrary distributions of v, respectively. In the case of the arbitrary distribution the 
value D 2 jV} = 0.9 has been chosen. 



Qt is shown in Fig. ||[ The curves g, p, and a correspond to the geometric, 
Poisson, and arbitrary distributions of u, respectively. In the case of the 
arbitrary distribution the value D 2 {z/} = 0.9 has been chosen. The sites of 
maxima of curves g, p and a are the same, but the heights of the maxima 
depend, however, on the character of the distribution of v. 

The relative variance of the number of living nodes at the time instant t 
can be written in the form: 



D 2 {M*)} 

t2 



(30) 



E ^>] q 2 Qt, if * 
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from which it is evident that the relative variance is increasing with t to 
oo in the subcritical and critical evolutions (qi < 1), and converging to the 
saturation value 

E{(^-1) 2 } 

if the evolution is supercritical (gi > 1). It is trivial but interesting to note 
that the standardized expectation value of fie(t) has the form: 



E{M*)> 
D{^(t)} 



i-gi 



1/2 r 



E{(^-1) 2 }, 



,[l-qi)Qt 



-1/2 



if 

if qi = l. 



(31) 



If the evolution is critical, then the standardized expectation value of the 
number of living nodes converges to zero as (Qt)^ 1 ^ 2 when t oo. 

4.1.2 Dead nodes 

It is obvious that the number of dead nodes Hd{t) in an evolving tree is a 



non- decreasing random function of time. From Eq. (|TTJ) we obtain for the 
expectation value E{/x<i(t)} = m,i(t) the following equation: 



mf\t) = 1 

the solution of that can be written in the form: 



,-Qt 



+ qiQ f V^-*') m[ d) (t')dt', 
Jo 



(32) 



mf\t) 



1-91 

Qt, 



--(l-qi)Qt 



if qi^h 
if qi = l. 



(33) 



We see that the expectation value of dead nodes converges to (1 — gi) -1 in 
subcritical and becomes infinite in critical and supercritical evolution when 
t =^ oo. It is remarkable that the increase of m^(t) with time t is linear in 
critical state. 

In order to calculate the variance D 2 {/i^(t)} of the number of dead nodes 
at the time instant t > we have to solve the integral equation 

m f{t) =q x Q C e -Q(*-*0 m ?(t') dt'+ 
Jo 
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Q jfe-<K*-*> \2q im f\t') + q2 [m ( *\t')} 2 
which can be derived from Eq. ( p]) by using the formula 



d 2 g^(t, 
dz 2 



z=l 



After elementary calculations we obtain the solution in the form: 

m f\t) = a { ? [1 - e~^ Qt }+ 
a? Qt e-^ Qt + e~^ Qt [1 - e^-^ Qt ], if Ql ^ 1, 



where 



2gi + 



92 



(l_ gi )2 i - !_ gi 



a? 



9i + 



92 



1 - gi V 1 - 9i 



and 



(d) _ 92 



If qi = 1, then we have 



(*) = (Q*) 2 + ~g 2 (Qt) 3 . 
By using the expression 

D 2 {Mt)} = rrtf>(t) + m?(t) [1 - mf\t)\ 
if qi ^ 1, then we find 

D 2 {M*)} = 4 d) [l-e-( 1 -^ t ] + 



(l-9i) 3 ' 



1 



4 d) Qt e^ 1 -^* + v { 3 d) e -^ Qt [1 - e~^ Qt ] 



where 



id) 



D 2 M , <0 _ „ D 2 M 



,(<*) 



E{(^-1) 2 } 



1 "(l-9i) 3 ' 2 (l-9i) 2 ' 3 (l"9i) 3 

and if gi = 1, i.e. the evolution is critical, then 

B 2 {fi d (t)} = Qt+^q 2 (Qtf. 
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One can see that the limit values of the variance of the number of dead nodes 
are given by the formula 



lim D 2 {/i d (t)} 

t — >oo 



4.1.3 Total number of nodes 



d 2 M :f „ ^ 1 




For many purposes, it may be important to know how the expectation value 
and the variance of the random function 



n(t) = nt{t) + n d (t) 



(40) 



depend on the time. By using Eq. ( |i5|) one can derive integral equations the 
factorial moments 



rrn{t) =E{p(f)} 
m 2 (t) = E{fi(tMt)-l}} 



dg(t, z) 
dz 

d 2 g(t, z) 
dz 2 



2 = 1 



2 = 1 



After a brief calculation we obtain that 

rt 

and 



mi (i) = 1 + qi Q f e -Q(*-*') mi (t') 
m 2 (t) = gi Q f e - Q W m 2 (f) df+ 



Q / ( 

/o 



mx(t) 



{2g 1 m 1 (t') + g 2 K(t')] 2 } df. 
The solution of Eq. ( fl3|) has the form: 

(1-gi)" 1 [l - qi e -^-^ Qt ] , ifgi^l, 

l + Qt if gi = 1, 

and looking at Eqs. (^) and ( |33|) we can see the trivial identity 

m 1 (t)=m[ e \t)+m{ d \t) 



(41) 
(42) 

(43) 



(44) 



(45) 
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to be valid. For the purpose of the calculation of the variance 



D 2 {/i(t)} = m 2 (t) + mx{t)[l - mi {t)} 



we have to obtain the solution of Eq. (HW). By having this solution we can 
immediately write: 

D 2 {/.(t)} = 

d 2 M 



d 2 M 



?1 



-(l-qi)Qt 



E{(^-1) 2 } 



(l-qi)Qt 



2qi 



1 



9i) : 



Qt e 



-(l-gi)Qt 



+ 



-(l-<7l)Qi 



if 91^1, 



(46) 



(l-9i) 3 

and in the case of critical evolution, i.e. if qi = 1, than we have a simple 
formula: 

D 2 {/i(t)} = (i + g 2 ) Qt + g 2 (Qt) 2 



1 



g 2 (Qty 



(47) 



The average tree size at a given time instant t > can be characterized by 
the expectation value of the total number of nodes /i(£), but it is to mention 
that the actual sizes of randomly evolving trees may fluctuate significantly 
around the expectation value E{yu(t)}. This fluctuation can be measured by 
the relative variance 

[EM()}] 2 

the time dependence of which shows interesting features. We see from Eq. 
(|45|) that the average tree size converges to infinity with t =>- oo , if q\ > 1, 



and to (1 — gi) , if (ft < 1. On the contrary, the relative variance remains 
finite when t =^ oo in the case of q% ^ 1, and if and only if the evolution is 
critical, i.e., if gi = 1, then can be observed the relative variance to become 
infinite when t =^ oo. 

In the part A of Fig. |] three curves of the relative variance of fi(t) versus 
Qt are plotted when qi = 0.95. In the part B the dependence of the relative 
variance at Qt = oo can be seen on the parameter q\. In both parts of Fig. 
the curves g, p, and a correspond to the geometric, Poisson, and the arbitrary 
distributions of u, respectively. At q± = 1 one can observe singularity in the 
relative variance, and thus one has to conclude that the average tree size at 
large Qt cannot characterize the random evolution in critical state. 
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Figure 4: Fig. A shows the relative variance of the total number of nodes vs. time 
parameter Qt. In Fig. B it can be seen at Qt = oo the dependence of relative variances on 
parameter q\. In both Figs, the curves g, p, and a correspond to the geometric, Poisson, 
and the arbitrary distributions of v, respectively. In the case of the arbitrary distribution 
D 2 -{V} = 0.9 has been chosen for the calculations. 



4.2 Covariances and correlations 

The study of how the random functions fie{t), fid{t) and /i e (t) are related 
involves the analysis of at least two correlation functions. In the following 
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we will investigate correlations between the living and dead nodes, as well as 
the lines and living nodes. 



4.2.1 Living and dead nodes 



In Sub-subsection 3.1.4 we derived an integral equation [see Eq. (|IBD] for 
the generating function g( e ' d \t, Z£, z d ) of the joint distribution p( e ' d \t, ri£, n d ) 
of random functions fJLi(t) and n d {t). It is clear that p^'^(t, ni, n d ) gives 
the probability that at time instant t > the number of living and that of 
dead nodes in a randomly evolving tree are equal to ri£ and n d , respectively, 
provided that at t = the tree was in the state So- 

Covariance function. In order to calculate the covariance function 



Cov{fx e {t)fx d {t)} = E{m{t)fi d (t)} - E{m(t)}E{fi d {t)} 
we need the equation determining the moment 

-d*g(W(t,z e ,z d y 



E{fi t (t)fi d (t)} 



J z e =z d =l 



dz e dz d 

The equation can be written in the form: 

mff{t) = Q qi f e ~ Q ^ m?f(t') dt'+ 



(48) 



(49) 



Q 



-Q(t-t') 



q x m?(t') + q 2 m ( p(t') (t')\ dt' . 



(50) 



By using the Eqs. (p6[) and fl33[ ) for rn(\t) and m}f\t), respectively, and 
introducing the notation a — (1 — qi) Q, it can be proved that 



Qt e~ at - q 2 



-<\t 



1 - e— ) 



is the unique solution of the integral equation (pOl) V qi ^ 1. After a brief 
calculation we have 



Cav{fM t (t)n d (t)} = Qt e" at - 



1 + 



D 2 M 

:i-<&) s 



1-e 



(51) 



if qi + I- 
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It is easy to show that in the case of critical evolution, i.e., if q\ = 1, then 

1 



(52) 



Remark Now, we would like to show that this result can be derived 
without solving the integral equation (|50|). Since = fii(t) + Hd(t) we can 
write that 

D 2 {/i(i)} = D 2 {^(t)} + D 2 {^(t)} + 2 Cov{^(t) fM d (t)}, 

where 

Cov{^(t) p d (t)} = E{[n £ (t) - m[ e \t)} \ji d (t) - m[ d \t)}}, 
and so we have 

Cov{^(t) = ~ [D 2 {/iW) - D 2 {/i,(t)} - D 2 {^(t)} 

By using Eqs. (BB), (E3) and (E|) we find 



Cov{ W (f) A*d(*)} = 0* e-d-^JO*. 

1 - gx 



D 2 {i/} " 

1 + 


e -(l-9l)Qt 


I _ e -(l-«i)Qt" 


L (i-?i) 2 . 







(53) 



if 

In the case of critical evolution by taking into account the expressions (^7), 
(130), and fl39T) the covariance can be written in the form: 



Cov{^(t) // rf (t)} = - g 2 (Qt) 2 . 



(54) 



These expression (|53|) is exactly the same as Eq. (|5l|) while (|54"D coincides 
with (0). 

The dependence of the covariance function 



Cav{» t {t) l i d {t)} = vV> i >{Qt, E{u}, D 2 {z/}) 
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qi 

Figure 5: Contour plot of the covariance function Cov{/^(t)/Aj(i)} = V^< d \Qt, qi) in 
the case of subcritical random evolution by assuming that the offspring number v is of 
geometric distribution. 

on the time parameter Qt and the expectation value as well as the variance of 
the v has been calculated by assuming the v to be of geometric distribution. 
Since in this case E{i/} = q x and D 2 {z/} — qx(l + qx), thus one can write 

The contour plot of this covariance function has been calculated in the 
(Qt,qi) plain for the subcritical values of q±. The results are seen in Fig. |5|. 
In order to show more precisely the form of time dependence of V^ l ' d \x, qi), 
three curves a, b, and c corresponding to the values q\ = 0.96, 0.97, and 
0.98, respectively, are plotted in Fig. |. One can observe that each of the 
covariance curves vs. time has a well defined maximum. The time parameter 
Qt m due to the maximum of the covariance function depends rather sensi- 
tively on q\. It is interesting to look at the variation of Qt m when q\ =>• 1. 
The curve plotted in Fig. [7| shows this variation, and one can see a steep 
increase of Qt m when q\ is approaching to 1. 

Correlation function. It is trivial that the numbers of living and dead 
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Figure 6: Three curves illustrating the time dependence of the covariance T>( e ' d ) (Qt, qi) in 
the case of subcritical random evolution by assuming that the v is of geometric distribution. 
The curves a, b, and c correspond to the values of q\ = 0.96, 0.97 and 0.98, respectively. 
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Figure 7: The dependence of Qt m on qi in the case of subcritical evolution assuming 
that the v is of geometric distribution. 

nodes are correlated, but the time dependence of the correlation function 

C (t ' gi ' 92) -D{Mt)}D{Mt)}' (55) 
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Figure 8: Time dependence of the correlation between the numbers living and dead 
nodes if v is of geometric distribution. The curves a, b, and c belong to q\ values 0.95, 1, 
and 1.05, respectively. 



has some specific features, namely, in the case of subcritical evolution the 
correlation function — after a sharp increase — is converging to zero when 
Qt =>- oo, while in the case of supercritical evolution reaching rapidly the 
value 1. 

If the random evolution is critical, i.e. q\ = 1, then we have 

Q 1-1/2 



C^(t,l,q 2 ) 



2 



1 + 



D 2 {z/} (Qtf 

which shows that the limit value of the correlation 

lim C^)(t,l,g 2 ) = ^ 



(56) 



(57) 



is independent of distributions of v and r. This very important limit law 
expresses the fact that the relation between the numbers of living and dead 
nodes near the critical state, and at time moments sufficiently far from the 
beginning of the process, becomes almost independent of rules controlling 
the evolution. 
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Figure 9: Dependence of the correlation between the numbers of living and dead nodes 
on qi at time instants Qt = 100, (a), Qtx = 500, (b) and Qt — oo, (c). 



Fig. |^ illustrates the time dependence of the correlation function 
, qi,q 2 ) in the case if v is of geometric distribution. The curves a, b, 
and c belong to the q\ values 0.95, 1, and 1.05, respectively. It can be shown 
that these curves are rather insensitive to the distribution of v. 

It seems to be worthwhile to present some curves reflecting the depen- 
dence of C^'^(t, qi, q 2 ) on qi at different time instants. In Fig. [] we can see 
that the transition of the correlation function from to 1 becomes sharper 
and sharper with increasing Qt. The curves a and b are belonging to the 
time parameters Qt = 100 and Qt = 500, respectively, while the curve c 
shows how the correlation function depends on q x at Qt = oo. One has to 
note that 



lim 

Qt->o 



C^ d \t,l + e,q 2 ) 



0. 



if e < 0, 



y/3/2, if e = 0, 



if e > 0. 
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4.2.2 Living nodes and lines 

In order to study the time dependence of correlation between the numbers 
of living nodes and lines in randomly evolving trees we should calculate the 
expectation value 

E{ yU ,(t) /Ue (t)} = m£ e) (t), 

which can be determined from the generating function g^' e \t, Z£, z e ) defined 
by (£?!]). Since 

,X dzg dz e 

one obtains from Eq. fl22|) the integral equation 

m [\ e \t) =Q qx f'e-Qit-t') m ff( t ') dt'+ 



- Z e =Z e = l 



Q I < 

10 



-Q(t-t') m (i)^ | gi + q2 \i + m [ e \t')] } dt'. 
Taking into account Eqs. 



-(l-qi)Qt 



and 



mfit) 



_^[l_ e -(i- 91 )Q*] ; if gi ^l, 



Qt, 



if qt = 1, 



it can be shown that the solution of (58) has the following form: 

9i 



i-9i (i - qi) 1 

if 9i ^ 1, 

where a = (1 — gi)<3, and 

m ^(t) = (l + g 2 )Qi + Ig 2 (Qt) 2 , if 



-at, 



1-e" 



91 = 1. 



By using the relation 

Cov{^(t) /i e (t)} = m£ e) (t) - m^t) m< e >(f) 



(58) 



(59) 



(60) 
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one can obtain that 



Cov{^(t) fi e (t)} = Qt e- a *- gi 



and 



1+ 1 { 



e -«*(l_ e -«t) > (6i) 



if qi 1, 



Cov{/i,(t) /i e (t)} = D 2 {z/} Qt (l + - Qi) , 
if 9i = 1. 



(62) 
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Figure 10: Time dependence of the correlation between the numbers of lines and living 
nodes in the case when v is of geometric distribution. The curves a, b, and c are belonging 
to the qi values 0.95, 1, and 1.05, respectively. 



Let us investigate the time dependence of the correlation function 



It can be proved that 



%(«)}=%(*)}, 



(63) 
(64) 
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and thus, if q\ ^ 1 then there is an easy task to calculate by using the 
expressions (|6T|), (|29|) and (|64|) , while if gi = 1 then we can write 



\/3 1 + 2 (Qt)- 1 

;i + 3(gt)-i(i + ^gfi(gt)- 1 



C l ; (t,l,g 2 ) = — 7 ; 77^7-, 777/2' ( 65 ) 



from that it follows the limit relation 



lim C^(t,l,q 2 ) 

Qt~*oo 2 

which is the same as lim^oo C™' d >(t, 1, q 2 ) and so, we can conclude that the 
limit value of the correlation between the numbers of lines and living nodes 
is completely independent of the parameters determining the tree evolution 
process. 

Fig. |lOj shows the time dependence of the correlation C^' e \t, qi, q 2 ) near 



the beginning of the process at the values of q± = 0.95, 1, 1.05 assuming that 
v has geometric distribution. The curve a corresponding to the subcritical 
evolution has a well defined maximum, and after that it decreases rapidly to 
zero when t =^> 00. In the case of the critical evolution, i.e., when qi = 1, the 
correlation function passing its maximum converges to the limit value v3/2 
as seen on the curve b. Finally, one can see that the curve c due to the 
value qi = 1.05 is reaching the level 1 not monotonously but through a slight 
maximum. The dependence of the correlation function C^' e \t,qi,q 2 ) on q± 
is similar to that of C^'^(t, gi, 52), i.e., the transition from to 1 becomes 
sharper and sharper with increasing Qt. 



5 Concluding remarks 

By introducing the notions of living and dead nodes a new model of random 
tree evolution with continuous time parameter has been constructed . In 
order to describe the evolution process two basic random functions [it{t) and 
Hd(t) have been used. The first is the actual number of living nodes, while 
the second that of the dead nodes at time instant t > 0. 

Having assumed the evolution process to be controlled by the lifetime 
and the offspring number of living nodes we derived exact equations for 
generating functions of fiiit), {^i(t), ^d(t)} etc., provided that at 
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t = the tree consisted of a single living node only. It is remarkable that 
the average lifetime has a role of scaling the time only in these equations. 

The time dependence of the average number of nodes has been used for 
characterization of the evolution which can be either subcritical or critical, or 
supercritical. It has been proved that the relative variance of the tree size is 
increasing with time, ie. at large t the average tree size is hardly informative. 
If the average offspring number converges to 1, then it has been found that 
the limit value of the relative variance tends to infinity. 

A specific property of the tree evolution has been discovered, namely, 
it has been proved that the correlation between the numbers of living and 
dead nodes decreases to zero in subcritical and increases to 1 in supercritical 
evolution, if the time parameter Qt tends to infinity, but in the case of 
exactly critical evolution it converges to a fixed value a/3/2 which is free of 
the parameters of the process. 
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